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ABSTRACT 

NGC4258 is the galaxy with the most accurate (maser-based) determination for the 
mass of the supermassive black hole (SMBH) in its nucleus. In this work we present 
a two-dimensional mapping of the stellar kinematics in the inner 3.0 x 3.0 arcsec 
= 100 X 100 pc of NGC 4258 using adaptative-optics observations obtained with the 
Near-Infrared Integral Field Spectrograph of the GEMINI North telescope at a « 0.11 
arcsec (4 pc) angular resolution. The observations resolve the radius of influence of the 
SMBH, revealing an abrupt increase in the stellar velocity dispersion within «10 pc 
from the nucleus, consistent with the presence of a SMBH there. Assuming that the 
galaxy nncleus is in a steady state and that the velocity dispersion ellipsoid is aligned 
with a cylindrical coordinate system, we constructed a Jeans anisotropic dynamical 
model to fit the observed kinematics distribution. Our dynamical model assumes that 
the galaxy has axial symmetry and is constrncted using the multi-ganssian expansion 
method to parametrize the observed surface brightness distribution. The Jeans dy¬ 
namical model has three free parameters: the mass of the central SMBH (M,), the 
mass-luminosity ratio (Tfc = M/L) of the galaxy and the anisotropy of the veloc¬ 
ity distribution. We test two types of models: one with constant velocity anisotropy, 
and another with variable anisotropy. The model that best reproduces the observed 
kinematics was obtained considering that the galaxy has radially varying anisotropy, 
being the best-fitting parameters with 3a significance M, = 4.8('(Q g x 10^ Mq and 
Tfc = 4 .I('(q' 5. This value for the mass of the SMBH is just 25 per cent larger than that 
of the maser determination and 50 per cent larger that a previons stellar dynamical 
determination obtained via Schwarzschild models for long-slit data that provides a 
SMBH mass 15 per cent lower than the maser value. 
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1 INTRODUCTION 

NGC 4258 (M 106) is a spiral galaxy w ith Hubble type 
SABbc, at a distance of 7.28±0.3 Mpc llHerrnstein et all 
Il999l) . which harbours one of the closest active galactic nu¬ 
cleus (AGN), classified as Seyfert 1.9. This galaxy is well 
known from previous studies for harbouring the supermas¬ 
sive black hole (SMBH) with the best constrained mass 
after that of the Milky Way, with a value of M. Maser = 
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(3.82 ± 0.2) X 10^ Mq, obtained from resolved kinematics of 
a rotating H 2 O maser disc within 0.13 pc from the nucleus 
l[Mj yoshi et al.iri995l : iGreenhill et al.|[l995l : iHerrnstein et al] 
I 1999 I) . It is also well known for its anomalous arms, which 
resemble spiral arms but are more di ffuse than regular spi¬ 
ral arms (IWilson, Yang CecillFoOll) . The anomalous arms 
span 5 arcmin in opti cal line emission a nd 10 arcmin in radio 
continuum emission (ICecil et al.ll2000l ). 

A nuclear non-stellar continuum and broad opti- 
cal emission l ines have been observed in polarized light 
llWilkes et al.l Il995h . supporting the presence of an ob- 
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scared AGN, from which a radio jet is observed propagat- 
ing p erpendicularly to the maser disc ll Cecil. Wilson fc Tullvl 
Il992h . The jet is oriented at a position angle (hereafter PA) 
PA=—3° ± 1°, consistent with the projected spin axis orien- 
tation of the mase r disc, which has a major axis PA=86°±2° 
(|Cecil et aLll2nnrit l. Farther from the nucleus, interactions be¬ 
tween the radio jet and the int erstellar gas are probably the 
origin of the anomalous arms dWilson, Yang &: Cecilll200ll . 
and references therein). 

In the infrared (hereafter IR), the nuclear continuum 
has been found to be well reproduced by a power-law /„ oc 
whic h seems to exten d through the optical to the 
ultraviolet (|Yuan et al.ll2003 l. In X-rays, the nuclear spec¬ 
trum presents several components: two power-l aws, a ther¬ 
mal component and the FeKa emission-line dYang et al.l 
l2007t l. The nuclear source is not resolved in the IR (1 to 
18/rm spectral region) in observations with 0.2 arcsec angu- 
lar resolution and presents total luminosity of 2 x 10® Lq 
l|Charv et aljlioool l. 

The host galaxy disc has the photometric major axis 
at PA =150°, as obtained from Hi observations llvan Albadal 
Il980ll , who have also obtained an in clination for t h e disc 
of i = 72°. A previous Ha study by Ivan der KruitI lll974| j 
gave values of PA = 146° and i = 64°. More recently 
ISawada-Satoh et al.l (l2007l l have mapped the large scale 
molecular CO(2-l) velocity field and concluded it is dom¬ 
inated by rotation in a disc with major axis PA=160° and 
1 = 65.6°. 

Other studies include the one of IPastorini et aO ll2007ll . 
who measured the gas kinematics close to the nucleus using 
long-slit spectra obtained with the Space Telescope Imag¬ 
ing Sp ectrograph (HST-STIS) and that of iHerrnstein et ^ 
1 I 2 OO 5 I I. who constrained even further the SMBH mass using 
ne w maser observatio ns. There is also a more recent study 
bv ISiopis et al.l 1 I 2 OO 9 II which uses the long-slit spectra from 
HST-STIS to obtain the stellar velocity distribution using 
the Can triplet absorption lines from 2 to 18.2 arcsec along 
the major axis and to 11.7 arcsec along the minor axis of the 
galaxy, obtaining a direct determination of the mass of the 
SMBH. 

Adopting the SMBH mass obtained from the kinemat¬ 
ics of the maser disc, and using a velocity dispersion for the 
galaxy bulge of lOSkms^'Q, one can calculate the radius 
of influence of the SMBH RilSpc. For a scale of 35 pc per 
arcsecond, this radius corresponds to 0.42 arcsec, thus re¬ 
solvable with ground-based adaptative optics observations, 
that typically reach 0.1 arcsec in the near-infrared. We thus 
decided to observe NGC 4258 with the Gemini Near-Infrared 
Integral Field Spectrograph (NIFS) in order to map the stel¬ 
lar kinematics and investigate the effect of the presence of 
the SMBH in the inner dynamics of the galaxy. Although 
we also use stellar kinematics t o pro be the presence of the 
SMBH as done bv ISiopis et al.l (I 2 OO 9 I I. our i^-band data are 
less affected by reddening, allowing the mapping of stellar 
kinematics much closer to the nucleus than previous observa¬ 
tions. In addition, our NIFS data provide a two-dimensional 


^ This is the average value of the velocity dispersion in the whole 
NIFS field of view and not the luminosity weighted velocity dis¬ 
persion inside a single aperture that contains the whole galaxy 
bulge (see Sec. IS. 
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Figure 1. Top Left: the NICMOS PSF modelled with the tiny 
TIM software. Top Right: the central 1.5 arcsec of the NICMOS 
image of the galaxy presenting the the features of the NICMOS 
PSF. Intensity units are normalized counts. The bottom panels 
show cuts along the two diagonals (dashed line and cyan contin¬ 
uum line) across the upper panels, in logarithmic scale. 


coverage instead of the limited spatial coverage of long-slit 
observations. 

Although the use of dynamical models based on the 
Schw arzschild orbit superposition method {| Schwar zschildl 
[l^ is currently the most popular method to determine the 
masses of SMBHs in the centre of galaxies, the Jeans dynam¬ 
ical models have proved capable of reproducing the values 
of the masses of the SMBH in good agree ment with those 
determined by the Schwarz schild models l|CaDDellari et all 
l2010l : lKrainovic et al.ll2009h . The Jeans models are also use¬ 
ful when one wants to have more physical insight to the 
problem. In addition, a Schwarzsc hild orbit superpo sition 
model of NGC 4258 already exists (|SioDis et akl 1200^ 1. pro¬ 
viding a value for the mass of the SMBH of M, schw = 
(3.3 ± 0.2) X IO^Mq, which is comparable to the better con¬ 
strained value (M.Maser = (3.82 ± 0.02) X lO’^M fT,) derived 
from maser observations IHerrnstein et al.l ll2005l l. 

We have thus decided to study the dynamics of the nu¬ 
clear region of NGC 4258 using our NIFS data and Jeans 
models in order to determine the dynamical quantities that 
govern the system. Our goal is also to investigate how well 
these models reproduce the observed stellar kinematics and 
how large are the differences between the values of the mass 
of the SMBH of NGC 4258 determined by different meth¬ 
ods. The use of the Jeans anisotropic dynamical model has 
allowed us to understand how variations in the model param¬ 
eters - such as the velocity anisotropy and the mass-to-light 
ratio - affect the stellar kinematics and the determination 
of the mass of the SMBH in NGC 4258. 

This paper is organized as follows. In section [2] we 
present our IFU observations and the photometric data used 
to construct the dynamical model. In Section [ 3 ] we discuss 
the derivation of the two-dimensional line of sight velocity 
distribution (LOSVD). SectionUis dedicated to the descrip¬ 
tion of the dynamical model and the discussion of the results 
and in section [5] we present our concluding remarks. 
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Figure 2. Top; the 2MASS ifs-band image. Bottom: the A-band 
image from NICMOS. The rectangle in the bottom panel shows 
the field of view of our NIFS observations presented in Fig|3] 


2 OBSERVATIONS AND DATA REDUCTION 

The Integral Field spectroscopic data were obtained in the 
near- IR A-band with the instrument NIFS ijMcGregor et ahl 
I 2 OO 3 II operating with the ALTAIR adaptive optics module 
at the 8-m Gemini North telescope in April 2007 under the 
science program GN-2007A-Q-25. The A-band was selected 
for this study because it contains the absorption CO bands 
at « 2.3 fim which can be used to derive the stellar kinemat¬ 
ics. In order to obtain the surface brightness distribution to 
model the gala:xy mass density distribution (and the gravi¬ 
tational potential) we used two sets of archive images: the 
first of them obtained with the Near-Infrared Camera and 
Multi Object Spectrometer (NICMOS) NIC2 camera of the 
HST. The second is a large scale image from the Two Micron 
All Sky Survey (2MASS). In the next sections we describe 
these sets of data and the reduction procedures. 


2.1 The photometric data 

2.1.1 NICMOS 

Near-infrared A-band images were obtained from the Mul¬ 
timission Archive at STScI (MAST) corresponding to data 
sets N46801050 for object observations and N46801060 for 
sky observations, from the proposal ID 7230 l|Scovi U3 ll997tl . 
These observations were obtained with the NIC2 camera 
of NICMOS using the F222M filter. The NIC2 camera has 
256 X 256 pixels, with scales of 0.075 arcsec pix“^ in x and 
y, providing a field of view of ~ 19.2 x 19.2 arcsec^. The 
F222M filter covers the wavelength range 2.15 — 2.28/rm, 
with central wavelength at 2.2174 fim. The galaxy observa¬ 
tion have a total exposure time of 224 seconds performed in 
MULTIACCUM mode with a four point spiral dither; the sky 
observations were obtained with a three point spiral dither. 

The data reduction was done with standard iRAij^ and 
PYRAF tasks following the same steps as in Chary et. al. 
(2000). We used the multidrizzle task to combine the 
dithered images following the standard sequence for NIC¬ 
MOS dithered observations described in the multidrizzle 
handbook. Inspecting the data quality files we noted that 
due the emission of the AGN, some of the pixels at the cen¬ 
tre of the galaxy are incorrectly identified as cosmic rays 
in the four dithered images. We thus corrected this in the 
data quality files and created a static mask to correct for a 
spurious “bar” feature in the column 128. Then a separated 
drizzled image was created for each dither position of the de¬ 
tector. To determine the need for extra offsets we have done 
cross correlation measurements between the images with the 
CROSSDRiz task and the extra shifts have them been com¬ 
puted with the SHIFTFIND task. Then we ran again the mul¬ 
tidrizzle with the extra shifts and created new separately 
drizzled images as well as a well-aligned median image. Fi¬ 
nally, we ran multidrizzle again to derive the cosmic ray 
masks and combine the images in a final drizzled image, 
clean of cosmic rays and detector artifacts. 

The compact nuclear emission due to the AGN its evi¬ 
dent in the reduced image. This unresolved emission appears 
in the image as a point source introducing the spurious fea¬ 
tures of the NICMOS point spread function (PSF) in the 
image. In the top left panel of Fig. [T] we show the NICMOS 
PSF, and in the right panel the central 1.5 arcsec of the im¬ 
age of NGC4258. In the bottom panels of this figure we 
display cuts along the two diagonals of each of the the im¬ 
ages in the top panels (black dashed and cyan continuous 
lines). The arrows indicate the strong signature of the PSF 
in the image of the galaxy. As our interest is to determine 
the luminosity distribution and gravitational potential gen¬ 
erated by the stars, we need to subtract the contribution 
from the AGN. In order to do this, we generated a large im¬ 
age of the NICMOS PSF and centred it at the same position 
of the galaxy centre in the detector. We then normalized the 
PSF image so that its peak flux coincided with that of the 
galaxy. After several tests in which we multiplied the nor¬ 
malized PSF by different fractions and subtracted it from 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 3. Left; Continuum map (reconstructed image) from our NIFS spectroscopy. Right; Spectra at the positions A, N and B marked 
at the left panel with the Bry emission-line and CO absorption bands identified. Flux units are ergs“^ cm“^ and the extraction apertures 
are 0.05 arcsec x 0.05 arcsec. 


the galaxy image, we concluded that the smoothest residual 
(without an “inverted peak” due to over-subtraction) was 
obtained with a contribution of the PSF of 70 per cent to 
the peak of the brightness distribution. The resulting im¬ 
age after reduction and AGN subtraction is shown in the 
bottom panel of Figure [2] The rectangle in the centre of 
the image represents the field of view of the NIFS obser¬ 
vations. The total amounth of light subracted corresponds 
to ~ 1.9 X10^ Lqk. In Apendix [C]we perform some tests to 
evaluate the effects of a possible over- or under-subtraction 
of the AGN on the determination of the mass of the SMBH. 


2.1.2 2MASS 

In order to characterize the surface brightness distribution 
of the galaxy beyond the inner 19.0 x 19.0 arcsec^ we used a 
large scale image fro m the 2MASS telesco pe in the A-short- 
band (As-band) from ljarrett et ^ (l2003l ). The 2MASS im¬ 
age has a scale of 1.0 arcsec pix“^ and the As-filter covers the 
wavelength range from 1.93 to 2.38/rm with central wave¬ 
length of 2.17 fim. The central 8.3 x 8.3 arcmin^ of this image 
is shown in the top panel of Fig. [21 The black rectangle in 
the centre of the panel represents the field of view of the 
NIGMOS image. 

2.2 The NIFS data 

NIFS has a square field of view of « 3.0 x 3.0arcsec^, di¬ 
vided into 29 slices with an angular sampling of 0.1 arcsec in 
X and 0.04 arcsec in y. The observing procedures followed the 
standard Object-Sky-Object-Sky-Object-Sky-Object dither 
sequence, with off-source sky positions since the target is 
extended, and individual exposure times of 600 s centred at 


A = 2.2 ym with a total of 10 individual exposures at the 
galaxy. The IFU was oriented with the slices along the po¬ 
sition angle PA= 145°. We have used the K_G5605 grating 
and the filter HK_G0603, which resulted in an arc lamp line 
full width at half maximum (FWHM) of 3.2 (thus R~5300). 

The data reduction was accomplished using tasks con¬ 
tained in the nifs package which is part of Gemini iraf 
package as well as generic iraf tasks. The reduction pro¬ 
cedure included trimming of the images, flat-fielding, sky 
subtraction, wavelength and s-distortion calibrations. We 
have also removed the telluric bands and flux calibrated 
the frames by interpolating a black body function to the 
spectrum of the telluric standard star. Small shifts between 
exposures due to guiding problems were corrected by mo- 
saicing the individual data cubes into a final one. The final 
data cube contains ~4300 spectra, each corresponding to an 
angular coverage of 0.5 x 0.5 arcsec^, which translates into 
1.75 X 1.75 pc^ at the galaxy. The NIFS angular resolution 
with ALTAIR in the A-band is 0.11 arcsec as measured from 
the FWHM of the spatial profile of the telluric standard 
stars. 

In the left panel of Fig.O we present a reconstructed 
image obtained from the data cube collapsing the spectra in 
the same wavelength range of the NIGMOS F222M filter, i.e. 
from 2.15 to 2.28/rm. The right panel of the figure presents 
three characteristic spectra extracted from the data cube: 
The nuclear spectrum (position N in the continuum map), 
a spectrum from a location at 1 arcsec NW of the nucleus 
(position A) and another from 1 arcsec SE of the nucleus 
(position B). The only emission-line present in the nuclear 
spectrum is the broad Bry at 2.1661/rm. The CO absorption 
bands used to obtain the stellar kinematics are identified in 
the spectrum from position B. 
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Figure 4. pPXF fits to the galaxy spectra. Examples of the re¬ 
sulting fit of the pPXf to the galaxy spectra on the positions A 
(top panel), N (middle panel) and B (bottom panel) indicated in 
Fig-El the cyan lines are the galaxy spectra, the black lines are 
the resulting fits of pPXF to the galaxy spectra and the red lines 
are the residuals. 


3 LINE OF SIGHT VELOCITY 
DISTRIBUTION 

In order to obtain the line-of-sight velocity distribution 
(LOSVD) w e have used the penalize d Pixel-Fitting (pPXF) 
methoc0 of ICappellari fc Emselleml ll2004h to fit the stel¬ 
lar absorptions bands present in the if-band spectra. The 
algorithm finds the best fit to a galaxy spectrum by convolv¬ 
ing template stellar spectra with the corresponding LOSVD. 
This procedure gives as output the radial velocity V, veloc¬ 
ity dispersion cr and higher-order Gauss-Hermite moments 
hs and h 4 . The pPXF method allows the use of several tem¬ 
plate stellar spectra and to vary the weights of the contri¬ 
bution of the different templates to obtain the best fit, min¬ 
imizing the template mismatch problem. We use the tem¬ 
plates from the spectroscopic library of late spectral type 
starfl observed with the Gemini Near Infr ared Spectrograph 
llWinge. Riffel &: Storchi-Bergmannl2009t) , which have an al¬ 
most identical spectral resolution to that of our data. We 
have used 60 templates in the pPXF fits, and for all spaxels 
of the Monte-Garlo simulations. We have also restricted the 


® Available from http://purl.org/cappellari/software 

Available from: http://www.gemini.edu/sciops/instruments/nearir- 
resources/?q=node/ 10167 


fit to the interval shown in Fig. 4, as for smaller wavelengths 
there are no features to constrain the stellar kinematics and 
there may be some contamination from emission lines. 


In Fig.|3] we present the resulting pPXF fits to the 
galaxy spectra from the positions A (top panel), N (mid¬ 
dle panel) and B (bottom panel) indicated in Fig. [31 The 
cyan lines are the galaxy spectra, the black lines are the 
pPXF fits to the galaxy spectra and the red lines are the 
residuals. The high signal-to-noise ratio of the IFU data, in 
the range of 30 to 130 with a average signal-to-noise ratio of 
~ 50, allows us to obtain fits of very high quality to the ob¬ 
served galaxy spectra without the need for spatial binning. 
The resulting kinematic maps of the velocity distribution 
are presented in Fig.O The top left panel shows the veloc¬ 
ity field after subtraction of the galaxy systemic velocity 
of 453kms“^, determined as the stellar velocity at the lo¬ 
cation of the peak of the continuum. The centroid velocity 
field presents a maximum of ~ 80 km s“^ and a rotation pat¬ 
tern with the line of nodes oriented along the position angle 
PA = 145° and with the SE side approaching and the NW 
side receding. The top right panel shows the map of the stel¬ 
lar velocity dispersion, varying in the range 80 — 180 kms“^, 
with a maximum at the location of the peak of the contin¬ 
uum. The average velocity dispersion over the whole field 
of view is ~ 105kms“^. In appendixlB] we make a compar¬ 
ison of t he resulting k inema tics with that presented in the 
paper of ISiopis et al.l l|2009l ). obtained with the HST/STIS 
instrument. As pointed out in the Introduction, the radius of 
influence of the SMBH in NGC4258 is Rinf ~ 15 pc, that cor¬ 
responds to 0.42arcsec at the galaxy distance. Fig.[3 shows 
that inside this radius there is indeed a sharp rise in the 
(T values as expected for the region where the SMBH dom¬ 
inates the gravitational potential, what is probably better 
seen in Fig.j^The lower panels show the Gauss-Hermite /13 
and /14 moments with values ranging from —0.15 to 0.15. 


The errors in the centroid velocities and velocity dis¬ 
persions were calculated from 850 Monte Garlo simulations 
and are illustrated in Fig. ( 6 ) The left panel of this figure 
presents the values of the measured centroid velocity (V) 
obtained from a cut along the galaxy major axis (black 
points) and the Icr error bar from the Monte Carlo sim¬ 
ulation (shaded band). The average standard deviation of 
the velocities is 3.3kms“^, with a maximum of 8.4kms“^ 
in the central pixel. Ignoring the region around the galaxy 
minor axis, where the velocities are approximately zero, the 
average velocity error is lower than 10 per cent. The right 
panel shows de values of the measured velocity dispersions 
(cr) (black points) and the Icr error bar (shaded band) along 
the galaxy major axis. The average standard deviation in 
the velocity dispersion is 4.1kms“^, with a maximum of 
10.2 per cent in the central pixel, an average error corre¬ 
sponding to 4 per cent. The errors in the central region of 
the data cube are higher than at the edges because of the 
presence of absorption features stronger in the nucleus than 
outside, that we attribute to the effect of dust and dilution 
by a red continuum from the AGN dusty torus. 
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Figure 5. Two-dimensional stellar kinematic maps of NGC4258. Top left; centroid velocity field. Top right: velocity dispersion. Bottom: 
Gauss-Hermite moments ha (left) and (right). Orientation as in Fig.[3] where the negative y-axis is at 145 degrees relative to the North 
direction. 


4 DYNAMICAL MODELLING 

4.1 Surface brightness distribution 

In order to model the galaxy surface brightne ss distribu¬ 
tion we use the mge_ fit-Sectors packag(0 (ICappellaril 
l2002t) which performs a Multi-Gaussian Expansion (MGE) 
parame trization for the observed surface brightness of the 
galaxy llEmsellem. Monnet fc Baconlll993) . In what follows, 
we present the most relevant steps in applying the MGE 
model to the galaxy NGC 4258 (see iGappellaril |2008l ) for a 
complete description of the method). 

The surface brightness distribution in the plane of the 
sky can be described by the sum of a set of two-dimensional 
concentric Gaussians as 


N 


k=l 




( 1 ) 


where N is the number of the Gaussians, each with total 


® Available from http://purl.org/cappellari/software 


luminosity Lk, axial ratio between 0 ^ g)) ^ 1 and disper¬ 
sion a'k along the major axis (a;'). In this reference system, 
the x' axis is aligned with the galaxy photometric major 
axis, derived from the photometry of the 2MASS image, 
which is oriented at PA=156°, and z' points to the line of 
sight. Before comparing © with the observed surface bright¬ 
ness distribution it is convolved with a MGE model for the 
NICMOS PSF. The NICMOS PSF w as obtained with the 
TINY TIM software (iKrist et al.ll201oh . The MGE parame¬ 
ters of the modelled NICMOS PSF are presented in Tab. lAll 
of ApoendixTAl In Fig. lAll we present a comparison of the 
MGE model and the tiny tim PSF. To examine how im¬ 
portant are the effects of the small differences between the 
MGE model and the tiny tim PSF in the convolution pro¬ 
cedure, we present in Fig. lA2l a comparison between the sur¬ 
face brightness distribution convolved with the tiny tim 
PSF and the one convolved with the MGE model for the 
PSF. Fortunately the two convolved surface brightness dis¬ 
tributions are very similar. This comparison show that the 
two are practically indistinguishable. 

In order to reproduce the surface brightness distribu¬ 
tion of NGC 4258 we fit together the 2MASS and NIC- 
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Figure 6. Measured values and Monte Carlo error estimation for the velocities and velocity dispersions along the galaxy major axis. 
Left: Measured values of the velocities V with the la error bar obtained from 850 Monte Carlo simulations. Right: Measured values of 
the velocity dispersions with the Icr error bar. 


Table 1. Gaussian parameters of the MCE model. 


I'j^ILqpc 2] 

[arcsec] 

<i'k 

247712.8 

0.063 

0.693 

28902.5 

0.280 

0.926 

9432.3 

0.501 

0.460 

10244.6 

0.981 

0.655 

3726.7 

2.677 

0.460 

2193.7 

3.525 

0.583 

607.5 

6.336 

0.460 

627.3 

9.022 

0.698 

387.3 

18.700 

0.658 

126.7 

41.920 

0.756 

181.2 

93.476 

0.450 

16.9 

247.722 

0.500 


Note: The first column lists the surface brightnesses of each Gaus¬ 
sian in units Lqpc“^; the second column list the Gaussian disper¬ 
sions along the major axis in arcseconds; the third column lists 
the axial ratios of the Gaussians. 

MOS A'-band images using the mge_fit_ sectors package. 
NGC 4258 has an angular size of 18.6 x 7.2arcmin^ . The 
NICMOS image has an angular size of ~ 19.2 x 19.2 arcsec^ 
and is used to fit the surface brightness distribution in the 
central region of the galaxy. Due to the large 2MASS PSF 
we do not use the the 2MASS image in the fit of the nu¬ 
clear region {R ^ 6.0arcsec). From the radius 6.0arcsec to 
the edges of the galaxy we fit the two images (NICMOS 
and 2MASS) together. As NGC 4258 is a spiral galaxy with 
a small bulge and a large disc component we implemented 
a two-component MCE parametrization by separating the 
Gaussians in two sets with different constraints for the ax¬ 
ial ratios {q'k). The first with 0.45 ^ q'k ^ 0.46 to model 
the disc component, and the second set with more flexible 
constraints, 0.5 ^ q'k ^ 0.99, to model the other galaxy 
components. 


After several attempts we obtained a satisfying fit of the 
galaxy surface brightness distribution using a set of 12 Gaus¬ 
sians centred in the galaxy nucleus and with the major axis 
aligned with the galaxy photometric major axis. In Tab.[T]we 
present the parameters of the Gaussians of our best fit: the 
first column lists the surface brightnesses (/(, = ) in 

units of L 0 pc“^, the second column lists the Gaussian dis¬ 
persions a'f. in arcsec and the third column lists the axial ra¬ 
tios q(,. In Fig.[7]we present linear cuts across the galaxy cen¬ 
tre: along the galaxy photometric major axis in the left panel 
and along the photometric minor axis in the right panel. 
The cyan continuous lines correspond to the MGE surface 
brightness distribution, the black open circles correspond to 
the observed brightness distribution in the NICMOS image 
and the red points correspond to the MGE after convolu¬ 
tion with the NICMOS PSF (E(g)PSF). In both panels the 
modelled convolved surface brightness distribution provides 
a good reproduction of the observed one for the inner 6 arc- 
seconds of the galaxy that is approximately two times the 
size of the region where we have kinematic measurements. 
Due the presence of the spiral arms and of isophotal twist 
the surface brightness distribution in the outer region of the 
galaxy can not be well reproduced by the MGE model. The 
good agreement of the red points with the open circles shows 
the good quality of the MGE model in describing the galaxy 
surface brightness distribution. 


4.1.1 Mass density and gravitational potential 

In order to obtain the deprojected three-dimensional lumi¬ 
nosity density we adopt the approximation that the galaxy is 
axisymmetric. In this case the luminosity density can be ob¬ 
tained from the parameters that describe the projected sur¬ 
face luminosity density. Assuming that the galaxy is oblate- 
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Figure 7. MGE model for the central 6.0 x O.Oarcsec^ of 
NGC4258 compared to the NICMOS image. Left: cuts along 
galaxy major axis of the modelled surface brightness distribu¬ 
tion (blue line), observed surface brightness distribution (black 
open circles) and the modelled convolved surface brightness dis¬ 
tribution (red points). Right: same as in the left for the galaxy 
minor axis. 


axisymmetric, the luminosity density can be obtained from 


V 






( 2 ) 


^q/2_j,Qg2 j 

where, qk = ^ ^ -, and i is the galaxy inclination 

(i = 90°being edge-on). It is important to make clear that 
this approach can not solve the intrinsic degeneracy of the 
deprojection. 

The galaxy mass density can be described by a set of 
Gaussians as 


must satisfy the steady-state collisionless Boltzmann equa¬ 
tion. Under the assumption that the galaxy has axial sym¬ 
metry we obtain t he axisymmetric Jeans equa tions in cylin¬ 
drical coordinates llBinnev &; Tremaind ll2008l l. equations 4- 
29). The solutions of these equations provide the velocity 
moments of the DF and can be compared with the observed 

kinematics of the galaxy. _ 

_ Semi-isotropic Jeans dynamical models djeand 


purposes ( 

Nagai & Mivamotol Il976l: ISatohl Il980l: 

Binney et al. 

I 990 I: van der Marel. Binney & Davies 1 

199d: 


Emsellem. Monnet & Bacon 1994: 

van Albada. Bertin & Stiavellil 1 19951: iRiciouti et al.l 

2 OO 5 II. Two 

of the most interesting applications are to 


estimate the d ynam i cal mass-to-light ratio of galaxies 
jvandei_^axel|_ 199ll;_|Stetlerj_Deion^he_&_Smecker;Han^ 


Cappellari et al.l 20061: Cortes. Kenney k. Hardy 


199C 


200811 and to obtain the masses of the SMBH in the nuclei 


of galaxies (iMagorrian et al.l 199^ Ivan der Marel et al.l 
1199^ CretJon_fc vjm den Bosch 19991 : 1 Joseph et al.l 2001^ . 

ICappellaril 1 2008l l has presented an anisotropic general¬ 
ization of the axisymmetric semi-isotropic Jeans formalism. 
Making the assumption that the velocity ellipsoid is aligned 
with a cylindrical coordinate system {R, z, (f)) and that the 
anisotropy is constant and quantified by v%, = bv'j, the Jea ns 
equations reduce to (eqs. (8) and (9) of ICappellaril ll2008tl l: 


)uv1-uvl d{bvvl) 


R 


dR 


'm 


(4) 


dz 


dz 


(5) 


p{R,z) 


M 


E 




(3) 


In the self-consistent case, i. e., if only stars contribute, we 
can obtain the galaxy mass density by the multiplication 
of the luminous density by a dynamical mass-to-light ratio. 
In this case the Gaussians in are the same as in © 
with Mj = Tfe Lfc. The gravitational potential generated by 
this density can be obtained through the Poison equation, 
= 47rGp. 

The contribution of the SMBH to the gravitational po¬ 
tential is modelled by the approximation that the corre¬ 
sponding distribution of matter is given by an extremely 
narrow spherical Gaussian. The dispersion of the Gaussian 
is constrained by the resolution of the kinematic data, such 
that 3cr. ^ Rmin, where Rmin is the smallest distance from 
the black hole that one needs to model. 


4.2 The Jeans anisotropic dynamical model 

A complet e description of the derivation of the Jeans 
equations lljean 2 1192211. can be found in the book of 
iBinnev &: TremaiM i 2008tl . Here we only provide an 
overview of some fundamental assumptions of the model. 
The model assumes that the galaxy is a large system of stars 
with positions (x) and velocities (u) described by a distribu¬ 
tion function (DF). If this system is in a steady state under 
the influence of a smooth gravitational potential, the DF 


One advantage of this approach is that, if a MGE model 
for the galaxy surface brightness is available, the solutions 
of the Jeans equations, (i. e., the first and second moments 
in velocity of the distribution function (DF) of the system, 
and the projections of these moments in the plane of the 
sky) are given as a function of the Gaussian pa rameters of 
the lu minosity density m and mass density Q . iGappell^ 
(l2008tl applied this method to determine the mass-to-light 
ratio and inclination of galaxies classified as fast-rotators 
in the SAURON survey Ide Zeeuw et al.ll2002l l. The effec¬ 
tiveness of this Jeans Anisotropic Multi- Gaussian expan¬ 
sion d ynamical model (JAM)0 is shown in iGappellari et ^ 
l|2009h where they model the observed kinematics of Gentau- 
rusA with the same parameters {M/L, M,) obtained from 
the Schwarzschild orbit superpositio n method. Another test 
of the JAM models is presented bv iLablanche et al.l (|2012h 
where the authors demonstrate the ability of the models 
in reproducing the anisotropy profile and mass-to-light ra¬ 
tios of realistic N-body co llisionless simu l ations of barred 
and unbarred galaxies. In iMedling et al.l 1 I 2 OIII I the JAM 
models are used to determine the mass of the SMBH in 
he southern c omponent of the pair o f inter acting galaxies in 
NGG 6240. In iNeumaver Sz Walcherl (I 2 OI 2 I') they established 
upper limits for the black hole masses of 9 l ate ty pe galax¬ 
ies using the JAM models. iR.aimundo et al.l ll2013l l uses the 
JAM models to obtain an upper limit of the black hole mass 


® Available from http://purl.org/cappellari/software 
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Figure 8. Determination of the best-fitting parameters and the minimization for the model C. Left: the plane M* x for b= 1.1 
(^f3z =0.09). Central: the plane M* x b for Tx =4.1. Right: the plane x b for M* = 4.7 x 10^ Mq. The values shown at the contours 
are for the quantity Ax^ = ~ Xmin where blue lines indicate the best-fitting parameters in each plane. The two 

inner contours shown correspond to the values ofAv^ = 2.3 and Ax^ = 11-8 that are the boundaries of the Icr and 3cr confidence regions 
for the two parameters together llPress et allJ ((20o3))- 


in MCG-6-30-15. lEmselleml (|2013l i used the JAM models to¬ 
gether with N-body simulations to review the hypothesis of 
the presence of an overmassive black hole in NGC 1277. An¬ 
other frequent application of the JAM models is the deter¬ 
mination of intermediate-mass black holes in globular clus- 


ole s m g l obular cius- 
aD l201ll.[2^.l2013l: 


ters llSeth et alJ 201^; iLiitzgendorf et aL ___ 

iFeldmeier et al.ll2013h . The largest applications of the JAM 
method t o date are the studies of the stellar initial mass 
function ijCap pellari et ^ 2012tl and dynamical scaling re¬ 
lations (|CaDDellari et al.ll2013l l of the 260 early-type galaxies 
of the ATLAS®° survey. 


The anisotropic Jeans method used in this work is sig- 
nihcantly different from Schwarzschild method. The former 
describes the orbital distribution via a few anisotropy pa¬ 
rameters, while the latter makes virtually no assumptions 
on the orbital distribution. A drawback of the Jeans ap¬ 
proach is that it can potentially bias the BH mass deter¬ 
mination. However this issue is largely overcome by allow¬ 
ing for a variation in the ansiotropy, which can be con¬ 
strained by integral-field kinematics, as we do here. Pre¬ 
vious tests agree with the result of this paper, showing a 
general consistency betw e en th e Schwarzsch il d and JAM ap¬ 
proach ICappellari et ah] ll201(]tl ; ISeth et ^ ll2014l l . An ad¬ 
vantage of the Jeans approach is its good predictive power: 
It is unable to fit the noise of the systematics in the data 
and can often be use d to detect prob l ems w ith the data 
or flag bad data (e.g. ICappellari et al.) (l2013lH . The Jeans 
solutions are also relatively easy to qualitatively under¬ 
stand. The Schwarzschild approach, given its generality in 
the adopted orbital distribution, does not suffer from po¬ 
tential bias in the BH mass. However the method can easily 
fit noise and bad data without raising any concern. The 
method can easily create unphysical orbital distributions 
ou tside the region constrained by the kinematics (e.g. Fig. 2 
of ICappellari &: McDerrnidI (l2005l H. For this reason it pro¬ 
vides better constraints to the SMBH masses when integral- 
field data are available out to larger radii (IVerolme et al.l 
I 2 OO 2 II . Overall, the two methods are complementary and 
sufficiently different to motivate the re-determination of BH 


mass presented in this paper, using Jean rather than the 
previously published Schwarzschild approach. 


4-2.1 The velocity second moment 

For a galaxy with a surface brightness distribution 
parametrized by a MGE model as in sec. 14.11 the projection 
in the plane of the sky of the second moment, , y')) of 

the DF is given in terms of t he parameters of the Gaussians 
(eq. (28) of ICappellaril ll2008h i. This quantity is a function of 
three free parameters: the galaxy mass-to-light ratio (Tk), 

the anisotropy {j3z = 1 ~ and the mass of the SMBH 

(M,). The comparison of the modelled with the 

measured Vrma = + a^, where V and a are shown in 

Fig-El provide the values of the best-fitting parameters for 
the galaxy. 

We start by assuming that the galaxy has a constant 
anisotropy in the velocity dispersions. Thus the space of pa¬ 
rameters has three dimensions: M, x Fk x /3z- We then per¬ 
form a minimization in this space searching for the values 
of M,, Fk and Pz that best fit the observed kinematics. We 
weight the minimization by assigning errors to the kine¬ 
matic measurements that are inversely proportional to the 
galaxy surface brightness at each location. This is done in an 
attempt to give comparable contributions to the x^ for all 
radii sampled by the NIFS kinematics. Without this weight¬ 
ing, the x^ would be artificially dominated by the numerous 
pixels at large radii, which contain virtually no information 
on the mass of the central supermassive black hole. In our 
trials without weighting the kinematics in the minimization 
we have obtained a significantly larger value for the mass 
of the SMBH (M. = 7.2 x 10^ M©), but the modelled val¬ 
ues of the velocity second moment clearly did not reproduce 
the values of the observed Vrms for the central region of the 
galaxy. 

We considered in a first trial that the galaxy has an in¬ 
clination of i = 72° (Model A). This inclination is obtained 
considering that the outer disc of the galaxy (r = 200 arcsec) 
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Table 2. Summary of the models with constant anisotropy. 


Model 

(1) 

i 

(2) 

M. |Mq] 

(3) 

Px 

(4) 

b m 

(5) 

■ 

'^min 

(6) 

A 

72° 

5.2x10'^ 

4.3 

0.90 (-0.11) 

27.21 

B 

68° 

5.0x10'^ 

4.3 

0.95 (-0.05) 

26.84 

C 

64° 

4.7x10'^ 

4.1 

1.1 (0.09) 

25.48 


Note. Column (1): model designation. Column (2): galaxy incli¬ 
nation. Column (3): best-fitting mass of the SMBH. Column (4): 
best-fitting /f-band mass-to-light ratio. Column (5): best-fitting 
anisotropy. Column (6): lowest value of the obtained. 

is thin (|van Albadal Il98fll : ISiopis et al.l l2009l l . The best¬ 
fitting model for this inclination with Xmin =27.21 is ob¬ 
tained with the parameters M, = 5.2 x 10^ Mq, rx=4.3 
and b = 0.90 (/3z =-0.11). This model does not reproduce 
satisfactorily the measured kinematics along the gala:xy mi¬ 
nor a:xis. 

After several attempts with different values for the 
gala:xy inclination we obtained a satisfactory model for i = 
64° (Model C) . The minimum obtained was Xmin = 25.48 
for M. = 4.7 X lO^'Me, rK=4.1 and b = 1.1 = 0.09). 

In Fig.| 8 ] we show the minimization in the space of 
parameters. In Tab. 2 we list the best-fitting parameters 
for three models with different inclinations: the models A 
and C described before and a third model with an inclina¬ 
tion of i = 68 ° for which the best-fitting parameters are 
M. = 5.0 X 10^ Mq, Tk =4.3 and b = 0.95 (/3^ =-0.05), that 
provide a minimum x^ of 26.84. Systematic variations in the 
best-fitting parameters are observed as the gala:xy inclina¬ 
tion decreases: the values of the mass of the SMBH and the 
mass-to-light ratio decrease but the value of the anisotropy 
parameters increase. 

In Fig.|^ we present the x^ minimization for M, for 
models A (black dotted line), B (open circles) and C (green 
solid line). The values of the parameters Fx and b are the 
best-fitting values of each model. The vertical line indicates 
the value of the maser determination, M, Maser = 3.82 x 
10^ M 0 . The green shaded region corresponds of the M, 
values that are within 3a conhdence intervals for the model 
with i = 64°. The best-fitting mass of the SMBH with 3cr of 
confidence is M, = 4.7lj g x 10^ M©. 

4-2.2 Effects of the PSF in the best-fitting parameters 

One important issue in the dynamical models of galaxies 
is that before comparing the modelled velocities with the 
measured ones we need to convolve the modelled veloci¬ 
ties with the PSF of the kinematic observations. For the 
Jeans anisotropic dynamical models the convol ution is im¬ 
plemented as in Appendix A of lCappellaril (l2008l l . We used a 
model for the NIFS PSF that is obtained from a MGE fit to 
the observations of three stars observed in the same night of 
the galaxy observations, as described in sec. IA2I As the time 
needed for the galaxy observations are longer (600 s) than 
that of the stars observations (15 s) there is the possibility 
that we are underestimating the NIFS PSF. 

In order to verify the influence of a possible broader PSF 
(larger FWHM) in the modelled kinematics and best-fitting 
value for M. we ran again the model C with Fx = 4.1, and 



Figure 9. Variation of xf' with the SMBH mass for the mod¬ 
els with constant anisotropy The green solid line represents the 
minimization of x^ as a function of M, for the best-fitting model 
with i = 64° (Model C). The values of M, in the shaded re¬ 
gion correspond to the 3 (t uncertainty in this parameter, thus the 
best-fitting mass of the SMBH is M, = 4.7^J g x 10^ Mq. The 
open circles represent the minimization for the model with i = 
68° (Model B). The black dashed line represents the minimization 
for the model with i = 72° (Model A). The vertical line indicates 
the value of the maser determination, M, Maser = 3.82 X 10^ Mq. 

b = 1.10 and performed the x^ minimization for M,. We 
used two broader PSF models obtained from the model de¬ 
scribed in sec. lA2l in the first one we increase the dispersions 
(listed in second column of Tab. lA2ll of the Gaussians of the 
MGE model of the PSF by a factor of 1.5 and in the second 
we duplicated the values of the dispersions of the Gaussians. 

In Fig.[T2]we present the modifications in the x^ mini¬ 
mization and in the best-fitting value of M, introduced by 
the increase of the NIFS PSF by a factor of 1.5 (represented 
by the black dotted line) and by a factor of 2.0 (represented 
by the open circles). The continuous line is for model C 
which provides a Xrnin = 25.48 for M. = 4.7 x 10^ Mq. 
For the model with the PSF 1.5 times broader we ob¬ 
tained Xmin = 26.17 for M. = 4.9 x 10^ Mq which is 
within the la uncertainties of the best-fitting model. For 
the model with the PSF 2 times broader Xmin = 33.83 for 
M, = 5.0 X 10^ Mq. For both cases the values for the mass 
of the SMBH that best reproduce the measured kinemat¬ 
ics are inside the 3a confidence interval of the best-fitting 
model and the modelled velocities still reproduce the ob¬ 
served ones. 
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Figure 11. Comparison between the observed Vrms and the best-fitting (Model D). Top left: the observed Vrms = \/V^~+~o^ 

after symmetrization; the contours correspond to the isovelocity curves of 104, 117 and 140kms~^. Top central: the best-fitting model 
presenting the superposition of the isovelocity curves of the modelled (red solid lines) and observed (black dotted lines) velocities. Top 
right: the residual Vrma — ™ kms“^. Bottom left: linear cuts along the galaxy major axis {x'), where the black points show the 

observed symmetrized Vrms, the red solid lines are the result of our best-fitting model, the shaded bands are the la errors in the velocity 
measurements. Bottom central and right panels show the same quantities as the left panel along the galaxy minor axis and along the 
diagonal direction, respectively. 


4-2.3 Models with variable anisotropy 


Table 3. Summary of the models with variable anisotropy 


As the solutions of the Jeans equations are presented for 
each Gaussian individually and the total second moment is 
the quadratic sum of the contributions of each Gaussian, it 
is possible to assign different values of anisotropy for each 
Gaussian of the MGE model that describes the ga l axy lu mi- 
nosity density. As exemplified bv iKrainovic et ^ (l200 9l) for 
the ga laxies NGG 2549 and NGC 524 and by Gebhardt et al.l 
l|200jh . spatial variations in the velocity anisotropy are fre¬ 
quent. In order to test how better the measured kinematics 
can be reproduced by models with variable anisotropy, we 
considered two scenarios. In the first we assumed that there 
is a radial variation in the anisotropy. We considered two dif¬ 
ferent values for the anisotropy parameters bk, one for the 
four inner Gaussians that have at < 1.0 arcsec and another 
value for the remaining outer Gaussians. Then we performed 
the minimization for two different inclinations i = 64° 
(model D) and i = 72° (model E). In a second scenario we 
assigned one anisotropy value for the Gaussians with ax¬ 
ial ratio Qk < 0.5 representative of the disc component and 
another anisotropy value for the Gaussians with axial ra¬ 
tio qk ^ 0.5 representative of the spheroidal component. 
This minimization was also performed for two values of the 
galaxy inclination: i = 64° (model E) and i = 72° (model G). 
Eor all the models with variable anisotropy we considered 
a constant mass-to-light ratio of Fk =4.1. The best-fitting 


Model 

i 

M. [Mq] 

bin Gz) 


■ 

■^min 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

D 

64° 

4.8x10’' 

1.10 (0.09) 

1.05 (0.05) 

25.24 

E 

72° 

5.3x10^ 

1.00 (0.00) 

0.90 (-0.11) 

28.65 




^disc 

^bulge i^z) 


F 

64° 

4.8x10’' 

1.25 (0.20) 

1.05 (0.05) 

25.34 

G 

72° 

5.3x10’' 

1.30 (0.23) 

0.85 (-0.17) 

27.90 

Note. 

Column 

(1): model 

designation. 

Column (2): 

galaxy 


inclination. Column (3): best-fitting mass of the SMBH. Column 
(4): for models (D) and (E) these are these are the values of 
the anisotropy parameter of the inner 4 gaussians, for models 
(F) and (G) these are the values of the anisotropy parameter of 
the Gaussians with axial ratio qk lower than 0.47. Column (5): 
for models (D) and (E) these are the values of the anisotropy 
parameter of the outer 8 Gaussians, for models (F) and (G) these 
are the values of the anisotropy parameter of the Gaussians with 
axial ratio qk higher than 0.5. Golumn (6): lowest value of the 
obtained. 


parameters for the models with variable anisotropy are pre¬ 
sented in Tab.O The best-fitting model with Xmin ~ 25.24 
was obtained for the inclination i = 64° and with a radially 
variable anisotropy (Model D). The dynamical parameters of 
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Figure 10. Effects of variations in the PSF width in the best¬ 
fitting mass of the SMBH. The green solid line represents the 
variation of with the mass of the SMBH for the model C. 
The dashed black line shows the variation of using a model 
for the NfFS PSF with a dispersion cr that is 50 per cent larger 
than that presented in App. m which results in a Xmin ~ 26.17 
for Mm = 4.9 X 10^ Mq. The open circles shows the variation of 
using a model for the NIFS PSF with a dispersion cr that is 
100 per cent larger than the original and results in a Xmin ~ 33.83 
for Mm = 5.0 X 10^ Mq. The shaded region shows the range of 
values of Mm that are inside the 3 cr uncertainty of the best-fitting 
value of Mm of the reference model (Model C). The vertical line 
indicates the value of the maser determination, Mm Maser = 3.82 X 
10'^ M0. 

this model are Mm = 4.8x 10^ Mq , Fk = 4.1, bin = 1.10 and 
bout = 1.05. The average absolute error in the second mo¬ 
ment over all pixels of the kinematic field is of 5.8 per cent. 
The results of the model for the second moment and the 
comparison with the observed Vrma are shown in Fig. llll 
The top left panel presents the observed Vrma after bi- 
symmetrization. The minimum value for Vrma is ~ 85 km s“^ 
and occurs at the most distant locations from the galaxy 
centre along the photometric minor axis. The maximum 
value is « 180kms“^ at the galaxy centre. The top cen¬ 
tral panel shows the distribution for the best-fitting 

model (model D), showing very good agreement with the 
measured Vrma- The top right panel shows the residuals 
Vrma — ■ The differences between the measured 

and modelled velocities are small, with the highest posi¬ 
tive residual being ~ 12kms“^ (which represents a error 
of « 11 percent in relation to the measured velocity in this 
position), and highest negative residual of « — 13kms“^ 



Figure 12. Variation of the x^ with the SMBH mass for the 
models with variable anisotropy. The green solid line represents 
the minimization of as a function of the Mm for the best-fitting 
model with i = 64° (Model D). The values of Mm in the shaded re¬ 
gion corresponds to the 3 ct confidence interval for this parameter, 
thus the best-fitting mass of the SMBH is M, = 4.8l!lg'9 ^ 

The black dashed line is the x^ minimization for the model with 
i = 72° (Model E). 


representing an error of « 12 per cent. The average error 
over the whole field is « 3 per cent. The good agreement 
between the modelled and measured velocities is also clear 
in the superposition of the isovelocity curves of the model 
on those of the measured velocity field presented in the top 
central panel. 

In the bottom panels we present linear cuts across the 
galaxy centre in three different directions: along the galaxy 
major axis in the left panel, the minor axis in the central 
panel and along a diagonal cut in the right panel. The black 
points show the measured values for Vrma with la error bars 
represented by the shaded region. The red solid lines show 
our best-fitting model for the projected second moment of 
the velocity {{vfoa)^^^)- The red dotted lines show the resid¬ 
uals between the observations and the best-fitting model. 
Along the galaxy major axis the modelled velocities repro¬ 
duced very well the observed ones. Along the galaxy minor 
axis the largest differences occur in the most distant regions 
from the galaxy centre. 

In Fig. [12] we present the minimization for M, for 
the models D (green solid line) and E (black dashed line). 
The vertical line indicates the value of the maser determi¬ 
nation, M. Maser = 3.82 X 10^ Mq. The green shaded re- 
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gion corresponds to the M, values that are inside the 3 (t 
confidence interval for the model with i = 64° (Model D). 
The best-fitting mass of the SMBH with Icr of confidence is 
M. = X IO^Mq. 

4 . 2.4 Comparison with previous results 

We can now compare the value of the SMBH mass that 
we have obtained with the previous ones from the lit- 
erature. When compared with the maser determination 
dCreenhill et aLlll99!a : iHerrnstein et al.l[i999ll of M, Maser — 
3.82 X 10^ M 0 our value of M, = 4.8lg g x 10^ Mq is « 
25 per cent larger. 

For comparison, the previous stellar dynamical deter- 
minatio n obtained via Sch warzschild models from long-slit 
data bv ISiopis et al.l ll2009l 'l [M, schw = 3.3±0.2 x 10^ M 0 ) is 
15 per cent lower than the maser value. Making a direct com¬ 
parison between the two dynamical determinations, there is 
a difference of « 45 per cent in the best-fitting values for the 
mass of the SMBH. The main factors that contribute to this 
difference are: 

• As we have shown in Fig. lBlI of appendixiBl the mea¬ 
sured velocity dispersions of our NIFS data are consider¬ 
ably higher than those tabulated for the STIS data, being 
the differences larger than the error bars for most of the 
points analysed. The differences between the values of the 
Vrms are consistent with the differences introduced in the 
modelled velocity second moment by the difference of the 
mass of the SMBH from the models. 

• Another possible cause for these difference in the mass 
of the SMBH can be an over-subtraction of the contribution 
from the AGN to the surface brightness distribution in the 
NICMOS image. The uncertainties introduced in the model 
by this factor are discussed in the appendixlcl 

• A third factor that could cause the difference in 
the SMBH mass derived using the two methods are 
the intrinsic differences between the dynamical mod¬ 
els used (Schwarszchild vs. Jeans models). The adopted 
JAM method makes more restrictive assumptions than 
Schwarszchild method. Alought the two methods have been 
shown to generally agree quite well in real galaxies, some 
differences are not unexpected. 

4 . 2.5 The veloeity first moment 

In order to model the velocity first moment it is necessary 
to make extra assumptions in the Jeans equations to specify 
how the second moment is composed in terms of ordered and 
random motions. We use the approximation that the tangen¬ 
tial component of the velocity first moment of each Gaussian 
is a function of the difference between the tangential and 
radial components of the velocity second moment of each 
Gaussian, ( yd,)k = K,k {{vfjk - (See Sec. 3.2.1 of 

ICappellm (l2008l l for a more complete explanation). 

Adopting the above assumption, we modelled the veloc¬ 
ity first moment using the best-fitting parameters of model 
E. A comparison between the observed rotation velocity field 
and the modelled first moment for this model is presented 
in Fig.lTSI 

The top left panel shows the measured velocity field for 
V. The top central panel shows the modelled velocity first 


moment (vios) with Kk = 0.92 and the best-fitting values of 
model D. The top right panel shows the residuals. The bot¬ 
tom right panel shows the residuals in percent. The bottom 
left and central panels compare the measured and modelled 
velocity fields along the galaxy major axis and along a di¬ 
agonal cut. It can be observed that the maximum residuals 
reach about 20kms“^ which correspond to « 50 per cent, 
as they are observed in the regions along the minor axis 
where the velocities are lowest. The residuals are higher 
than 30 percent within the inner 0.5 arcsec. In the outer re¬ 
gions the differences between the model and measurements 
are lower than 30 per cent. But the bottom panels highlight 
that the modelled velocity first moment does not reproduce 
the observed rotation mostly in the region between 0.1 and 
0.8arcsec (3.5pc and 30pc). 

The discrepancy between the modelled velocity first mo¬ 
ment field and the measured centroid velocity field, does not 
affect our conclusion about the mass of the SMBH, that is 
mostly based on the fit of the second moment. The first mo¬ 
ment derived from the model depends on an extra assump¬ 
tion: that the tangential anisotropy is the same everywhere, 
what may not be true. A more likely explanation for the 
discrepancy is however the inaccuracy in the deprojection 
of the luminous stellar density from the observed surface 
brightness. The deprojection is known to be mathemat i- 
cally degenerate dRvbickil Il987l : iGerhard Binti^ Il996ll , 
with the effect becoming quite sig nificant at low inclination 
llRomanowskv &: KochaneH Il997l l . Our galaxy is not seen 
close to edge-on and could easily hide a weak nuclear disc 
which would be invisible in the projected photometry, could 
be showing up in the velocity field, leading to the observed 
discrepancy between the model and the data. This state¬ 
ment is based on numerical experiments where we tried to 
force the model to have a flat disk by making the Gaus- 
sians of the MGE to have the same axial ratio. Only in this 
way one can lower the model inclination and make the in¬ 
ner disk intrinsically very thin within the NIFS FoV. With 
this forced MGE one can better reproduce the fast increase 
in rotation and the difference in Vrms between the major 
and minor axis. Unfortunately the galaxy overall is not well 
reproduced by an MGE with constant projected ellipticity. 
For this reason our test only shows that the existence of a 
very thin disk within the NIFS FoV would improve the fit 
to the kinematics. 


5 SUMMARY AND CONCLUSIONS 

We have presented a two-dimensional mapping of the stellar 
kinematics of the inner 3.0 x 3.0 arcsec^ (~ 100 x 100 pc^) 
of the galaxy NGC4258 using NIFS A-band data with a 
signal-to-noise ratio higher than 50 over most of the ob¬ 
served field, spectral resolution of 5300 and spatial resolu¬ 
tion of «4pc, allowing to resolve the radius of influence of 
the SMBH (« 15pc). The centroid velocity field presents a 
rotational pattern with a maximum velocity of ±80kms“^, 
with the SE side approaching and the NW side receding. The 
stellar velocity dispersion presents an abrupt increase within 
the inner 0.3 arcsec (10 pc), reaching a value of 180 km s“^ at 
the nucleus, consistent with the presence of a SMBH there. 

In order to model the stellar kinematics we built a 
Multi-Gaussian Expansion Jeans Anisotropic Dynamical 
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Figure 13. Comparison between the measured centroid velocity V and the first moment (vi^s) for model D. Top left: measured velocity 
V', the contours are for the isovelocity curves of 0.0, ±20, ±40, it60kms“^. Top central: modelled (ujog) with the best-fitting values of 
Mm = 4.8 X 10^ Mq, Tjy = 4.1, bin = 1-10 and bout = 1.05 obtained from the fit of the second moment and ft = 0.9. Top right: residuals 
in kms“^. Bottom right: errors in percent. Bottom left: linear cuts across the galaxy centre along the galaxy major axis. Bottom central: 
linear cuts across the galaxy centre in a diagonal direction. 


Model. In this model the velocity second moment {vfos)^^^ is 
a function of three free parameters: the galaxy mass-to-light 

ratio Fk, the anisotropy in the velocity b — = and the 

vi 

mass of the SMBH M,. We performed a minimization 
in this space of parameters in order to search for the values 
that best reproduce the measured Vrma = + a^). In 

order to ensure that all regions of the galaxy, independently 
of its radial position, have approximately the same relevance 
in the minimization process we have weighted the x^ mini¬ 
mization by assigning errors to the kinematic measurements 
that are inversely proportional to the galaxy surface bright¬ 
ness at each location. Without this weighting the best-fitting 
model does not reproduce the measured Vj-ms in the nuclear 
region of the galaxy and provides a a wrong (too high) value 
for the mass of the SMBH. 

We have tried models with only one value for the ve¬ 
locity anisotropy and models with two values. The best¬ 
fitting model was obtained adopting a galaxy inclination 
i = 64° and considering that the galaxy has a radially vari¬ 
able anisotropy in the velocity second moment, being the 
values of the anisotropy parameter for the inner 4 Gaus- 
sians bin ~ 1.10 {Pz ~ 0.09) and for the remaining 8 outer 
Gaussians bout = 1-05 {Pz = 0.05). 

Gonsidering the 3cr confidence intervals, we have ob¬ 
tained a mass-to-light ratio of F*, = d.l^Q g and the best¬ 
fitting SMBH mass of M, = d.S^Q g x 10^ Mq. These 3a 
uncertainties of our model are comparable to the typical 


values for the uncertainties of stellar dynamical models for 
other galaxies, that usually are lower than 50 per cent for the 
mass of the SMBH iMcConnell et all 1 I 2 OIIII . This provides 
further confirmation of the robustness of the stellar dynam¬ 
ical determinations. It is worth noting that an accurate BH 
masses was obtained even when using a simple model and 
only fitting the kinematics within the innermost few arcsec- 
onds, without the need for a large-scale model of the galaxy 
dynamics. 
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Table Al. Gaussian parameters of the MGE model for the NIG- 
MOS PSF. 


[Total Counts] 

CT^ [arcsec\ 

-Ifc 

0.52 

0.084 

Ri 1.0 

0.34 

0.242 

Ri 1.0 

0.06 

0.789 

Ri 1.0 

0.08 

2.671 

Ri 1.0 



Figure Al. NICMOS PSF fit. Left: linear cut along the detector 
horizontal direction showing the radial profile of the NICMOS 
PSF of the TINY TIM model (black open circles) and of the result¬ 
ing MGE model for the PSF (blue solid line). Right: Same as the 
left but for the detector vertical direction. 

APPENDIX A: MGE PARAMETERS FOR THE 
NICMOS AND NIFS PSF 

In this appendix we present the MGE models for the NIC¬ 
MOS and NIFS PSF’s and the agreement of the MGE model 
of the NICMOS PSF with the on of the tiny tim model. 

Al The NICMOS PSF 

We modelled the NICMOS PSF by a set of four approxi¬ 
mately circular Gaussians using the MGE method. The re¬ 
sulting Gaussian parameters are presented in Tab. \M\ In 
Fig.im we present the comparison between the luminosity 
profiles of the tiny tim PSF showed by the black open cir¬ 
cles and our MGE modelled PSF showed by the blue solid 
line. 

In order to verify the effects of the small differences be¬ 
tween the MGE model and the real PSF in the convolution 
procedure we present in Fig. lA2l a comparison of the con¬ 
volution of the surface brightness distribution of NGC4252 
presented in sec. 14.II with the tiny tim PSF (black open 
circles) and with our MGE model for the NICMOS PSF 
(red solid line) for the central 6.0 arcsec of the galaxy. The 
left panel shows a linear cut across the galaxy centre and 
along the photometric major axis and the right panel along 
the galaxy minor axis. In both directions the differences are 
very small in most of the regions. 

A2 The NIFS PSF 

The MGE model of the NIFS PSF is performed using the 
average image of three reconstructed images from the data 
cubes of three stars observed with NIFS in the same night of 



Figure A2. Comparison of the convolution of the MGE model 
for the surface brightness distribution with the tiny tim PSF and 
with the MGE model for the PSF. Left: cuts along galaxy major 
axis. Right: cuts along galaxy minor axis. 


Table A2. Gaussian parameters of the MGE model for the NIFS 
PSF. 


[Total Counts] 

[arcsec] 

I'k 

0.29 

0.045 

Ri 1.0 

0.20 

0.069 

Ri 1.0 

0.15 

0.103 

Ri 1.0 

0.27 

0.317 

Ri 1.0 

0.09 

0.543 

Ri 1.0 


the observations of the galaxy. The Gaussian parameters of 
the MGE model for the NIFS PSF are presented in Tab. lA2l 
In Fig. lA3l we present a comparison between the resulting 
MGE model and the NIFS PSF. In the left panel we show a 
linear cut along the detector horizontal axis of the NIFS PSF 
(black open circles) and of the MGE model (blue solid line). 
In the right panel we show a linear cut along the detector 
vertical axis. 



Figure A3. NIFS PSF fit. Left: linear cut along the detector 
horizontal direction showing the radial profile of the NIFS PSF 
obtained from observations of three stars (black open circles) and 
of the resulting MGE model for the PSF (blue solid line). Right: 
Same as the left but for the detector vertical direction. 
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Figure Bl. Comparison between STIS and NIFS kinematics 
along the galaxy major axis. Top: The blue diamonds are the 
measures of the centroid velocities from our NIFS data with the 
Icr error bar represented by the shaded region, the black circles 
with error bars are the measurements from STIS observations. 
Middle: Velocity dispersions from the NIFS and STIS observa¬ 
tions. Bottom: The blue diamonds and black circles are the NIFS 
and STIS measured Vrms respectively. The dashed line is our 
best-fitting model for the velocity second moment, the continu¬ 
ous line is the modelled velocity second moment using the mass of 
the SMBH from the maser determination and the d ashed-dotted 
line is the model using the mass of the SMBH from lsiopis et ahl 
ll2009h . 


APPENDIX B: COMPARISON BETWEEN THE 
KINEMATICS MEASUREMENTS OF NIFS 
AND STIS 


In this section we present a comparison between the 
measured velocities and velocity dispersions of our NIFS 
data and the previous publis hed long-slit data from the 
HST/STIS instrument used bv ISioois et alJ (l200flh to deter¬ 
mine the mass of the SMBH in the nucleus of NGC4258 
through Schwarzschild dynamical orbit superposition. In 
Tab.lBl] we present kine matic measurement s with la er¬ 
rors from the STIS data l|Siopis et al.1 l2009ll for the posi¬ 
tions along the galaxy major axis coincident with our NIFS 
measurements. In column one we list the positions of the ex¬ 
tractions; in column two, the values of the centroid velocities 
obtained from the STIS data; in column three, the values of 
the centroid velocities from the NIFS data; in columns four 
and five we present the values of the velocity dispersion mea¬ 
surements from the STIS and NIFS data respectively. 

The differences between the measured velocities from 
both instruments are shown in Fig. lBlI The top panel shows 
the centroid velocities: the blue diamonds are the NIFS data 
with the la errors being represent by the shaded grey band, 
the open circles with the error bars are the STIS data. In the 
region inside 0.2 arc seconds the STIS data shows a steeper 
velocity curve than the NIFS data, while for the outer re¬ 
gion the error bars overlap. We attribute this difference to 
the somewhat better angular resolution of the STIS when 
compared to that of NIFS. 

The middle panel shows the velocity dispersions: the 
NIFS data have cr values «20 km s“^ higher than 
those from the STIS data. We attribute this difference 
to the smaller FWHM of the STIS PSF and to the dif¬ 
ferent methods and templates used for the cr determina- 
tio ns: while we have used the stell a r tem pl ates library 
of Winge. Riffel &: Storchi-Bergmamil (120091 '). ISiopis et al.l 
(|2009ll used only one stellar template (HR 7615). 

In the bottom panel we show the resulting Vrms com¬ 
pared to different models. The black dctshed line is the ve¬ 
locity second moment along the galaxy major axis from our 
best-fitting model with the parameters bin = 1.10 = 

0.09), bnut = 1.05 (,5^ = 0.05), Fk = 4.1 and M. = 
4.8 X 10^ Mq; the black continuous line is the second mo¬ 
ment from a model with the same anisotropy and mass-to- 
light ratio of the best-fitting model but using the value of 
M, Maser = 3.82 X 10^ Mq ; and the red dashed-dotted line is a 
model using the value of M, schw = 3.3 x 10^ M©. While our 
best-fitting model gives a M, somewhat larger t han that of 
the maser determination, the lSiopis et all ll2009lf gives a M, 
somewhat smaller, but the NIFS measurements gives Vrms 
values closer to the one from the model with M. = M.Maser. 


APPENDIX C: SUBTRACTION OF THE 
POINT SOURCE 

In section ETTI we have determined that the luminosity of the 
AGN contributes with approximately 70 per cent to the total 
intensity of the central pixel of the NICMOS image. Thus the 
total luminosity subtracted from the galaxy image in order 
to eliminate the contribution of the AGN is of the order of 
~ 1.9 X10^ Lqk. Using the mass-to-light ratio of Fj, = 4.1 as 
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Table Bl. Comparison between NIFS and STIS kinematic measurements along the galaxy major axis. 


x' [arcsec] 

VsTlS [km s 1] 

VNiFs[fcms 1] 

o-STis[fc"is 

O-NIFS [ferns 1] 

0.00 

20± 14 

0±8 

159 ± 17 

183 ± 10 

0.05 

42 ±6 

17 ±8 

142 ±8 

170 ± 10 

0.10 

60 ±4 

31 ±6 

133 ±5 

155 ±7 

0.17 

59 ±4 

40 ±3 

122 ±5 

152 ±3 

0.30 

56 ±6 

51 ±2 

112 ±7 

120 ±3 

0.50 

57 ±7 

55 ±3 

85 ±6 

115 ±3 

0.80 

53 ±8 

62 ±3 

96 ± 10 

104 ±3 


Note: Column 1 are the positions o f the extract i ons al ong the galaxy major axis; Column 2 are the centroid velocities from STIS 
observations as tabulate in Table 3 of lSioois et al.l i|20Q9lj : Column 3 are the centroid velocities from the NIFS observations; Column 4 
are the velocity dispersions from the STIS observations; Column 5 are the velocity dispersions from the NIFS observations. 




Figure Cl. Effects of the point source subtraction on the modelled second moment - Test 1. Left panel: The continuous line is a cut 
along the galaxy major axis of our best-fitting model, the dotted line shows the modelled second moment using the the same dynamical 
parameters of our best-fitting model and subtracting 10 per cent (Ri 0.27 X 10^ Lq) more luminosity in the AGN contribution. The black 
fil circles represent the resulting model subtracting 10 per cent less luminosity in the AGN subtraction. In the two cases the luminosity 
is subtract from the innermost gaussian of the MGE model. The shaded band are the la confidence intervals for the Yrms- Right panel: 
Same as the left panel but for the galaxy minor axis. 


obtained from our best-fitting model this corresponds to a 
mass of « 7.8 x 10^ Mq. In order of investigate how a possible 
over- or under-subtraction of the AGN contribution affects 
the determination of the mass of the SMBH we performed 
two simplified tests. 

In the first test we assume that only the innermost 
gaussian component of the MGE model is affected by the 
subtraction of the point source. This is a plausible assump¬ 
tion as all the other Gaussians components of the model 
extend beyond the first Airy ring of the PSF and after the 
convolution of the surface brightness distribution with the 
PSF the Gaussian components become even flatter. We then 
added and subtracted « 0.27 x 10^ Lqk (or « 1.12 x 10^ Mq 
) from the innermost gaussian. In practice this is equivalent 
to instead of subtracting 70 per cent of the intensity of the 
central pixel to subtract 60 per cent and 80 per cent respec¬ 
tively. The effects on the modelled velocity second moment 
for the three different subtractions of the AGN contribution 
are shown in Fig. icn that presents cuts along the galaxy 
major and minor axis showing the variation of the modelled 
velocity second moments. The continuous line is our original 


best-fitting model (Model D), the dashed line is a model with 
the same dynamical parameters of the best-fitting model 
but with a additional mass of « 1.12 x 10^ Mq to the central 
gaussian, and the shaded band is the Icr confidence region for 
the kinematic measurements. The black filled circles are the 
model where we subtracted the same mass from the inner¬ 
most gaussian. The resulting variation in the velocity second 
moment of the central pixel of the galaxy for this variation 
in the mass is of the order of ±12kms“^. 

In the second test we assumed that the innermost gaus¬ 
sian has a luminosity of « 0.27 x 10^ Lqk higher, (this cor¬ 
responds to subtract only 60 percent from the intensity of 
central pixel, this test is motivated by the fact that this 
amounth of light corresponds approximately to a stellar 
mass that is equal to the difference in the values for the mass 
of the SMBH obtained from our best-fitting model and the 
value from the maser determination) and we use the param¬ 
eters of anisotropy and mass-to-light ratio of our best-fitting 
models to investigate how well the three different determi¬ 
nations for the mass of the SMBH (Maser, STIS, ours) re¬ 
produce the observed kinematics. The resulting modelled 
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Figure C2. Effects of the point source subtraction on the modelled second moment - Test 2. Left; the continuous line is a cut along 
the galaxy major axis showing the values of the velocity second moment for the parameters of our best-fitting model and subtracting 10 
per cent less luminosity (RiO.27 X 10^ Lq ) in the AGN contribution, The dotted line corresponds to the model using the value of the 
mass of the SMBH of the maser determination (3.82 X 10^ Mpi) and the black bullets is the model using the value of the mass of the 
SMBH from the previous stellar dynamical determination of lSioois et al.l (l2003h (3.3 X 10^ M 0 ).The shaded band are the Icr confidence 
intervals for the Vrms. 


velocity second moments along the galaxy major and minor 
axiz are presented respectively in the left and right panels of 
Fig. lC2l The continuous line corresponds to our best-fitting 
model, the dashed line is the model using the value for the 
mass of the SMBH from the maser determination and the 
black filled circles are the m odel with the value for the mass 
of the SMBH determined bv ISionis et al.l (l200fll '). The shaded 
band is the Icr confidence region for the kinematic measure¬ 
ments. Under these assumptions, the best-fit is obtained for 
the model with M.Maser = 3.82 x 10^ M©. The model using 
a mass of M, stis = 3.3 X 10^ Mq for the SMBH shows a 
poor but a similar fit to the data to that of our model only 
that our model a re close to the upp er envelope of the Icr 
band while the of lSioois et al.l fepod l are close to the lower 
envelope. 

But we would to point out that the uncertainty in the 
determination of the AGN is snaller than 10 per cent. Actu- 
aly even an over-subtraction of « 5 per cent already leaves a 
signature on the NICMOS PSF in the image. Our choice of 
« 5 per cent is because it corresponds im mass to the differ¬ 
ence between the value of the SMBH mass of our best-fitting 
model and that of the maser determination. 

































